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Abstract 

NASA has teamed with the FAA, DoD, industry, and academia for research into the remote detection 
and measurement of atmospheric conditions leading to aircraft icing hazards. The ultimate goal of this 
effort is to provide pilots, controllers, and dispatchers sufficient information to allow aircraft to avoid or 
minimize their exposure to the hazards of in-flight icing. Since the hazard of in-flight icing is the outcome 
of aircraft flight through clouds containing supercooled liquid water and strongly influenced by the 
aircraft’s speed and configuration and by the length of exposure, the hazard cannot be directly detected, 
but must be inferred based upon the measurement of conducive atmospheric conditions. Therefore, icing 
hazard detection is accomplished through the detection and measurement of liquid water in regions of 
measured sub-freezing air temperatures. The icing environment is currently remotely measured from the 
ground with a system fusing radar, lidar, and multifrequency microwave radiometer sensors. Based upon 
expected ice accretion severity for the measured environment, a resultant aircraft hazard is then 
calculated. Because of the power, size, weight, and view angle constraints of airborne platforms, the 
current ground-based solution is not applicable for flight. Two current airborne concepts are based upon 
the use of either multifrequency radiometers or multifrequency radar. Both ground-based and airborne 
solutions are required for the future since ground-based systems can provide hazard detection for all 
aircraft in airport terminal regions while airborne systems will be needed to provide equipped aircraft 
with flight path coverage between terminal regions. 

1.0 Introduction 

Aircraft Icing is a significant flight hazard that continues to be a major concern to the flying 
community (Ref. 1). A recent search of the Aircraft Owners and Pilots Association (AOPA) Accident 
Analysis Web site (Ref. 2) produced 62 structural in-flight icing accident reports covering the period from 
January 1999 to December 2007, where about half of these accidents resulted in fatalities. The severity of 
icing conditions can vary considerably both spatially and temporally, so that one aircraft may encounter 
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no perceptible icing while a trailing aircraft may encounter severe icing (Ref. 3). These factors conspire to 
make the accurate prediction of in-flight icing conditions very important, yet very difficult. While great 
progress has been achieved recently in the diagnosis and prediction of icing based on meteorological 
models and station reports (Ref. 4), the demand for the direct detection and measurement of hazardous 
flight icing conditions continues. 

Aircraft ice accretion typically occurs when an aircraft surface impacts supercooled liquid water. 
Depending on the amount of liquid water present, the size of the cloud droplets, the air temperature, the 
aircraft surface geometry, and the aircraft speed, the water will freeze immediately (typically lower water 
contents, smaller droplet sizes, and colder temperatures) which results in a generally conformal rime ice, 
or will flow on the surface before freezing (typically higher water contents, larger droplets, and warmer 
freezing temperatures) which results in the more hazardous, nonconformal glaze ice. Since aircraft icing 
is a process that requires the presence of an aircraft to occur, it cannot be directly measured remotely. 
Instead, it is the detection and measurement of the conditions that lead to aircraft icing in which we are 
interested. So, when we speak of the remote detection of aircraft icing, we are really speaking of the 
detection and measurement of supercooled liquid water. The current techniques for the remote detection 
and measurement of icing conditions rely on the identification and measurement of liquid water and then 
the inference or measurement of the surrounding air’s temperature. 

To provide coverage for the maximum number of aircraft over the maximum extent of their flight 
path, both ground-based and airborne icing remote sensing solutions are desired. While ideally all aircraft 
would be equipped with icing remote detection equipment to provide information coverage over their 
entire flight paths, this will practically never happen. So providing a ground-based capability that can 
provide information to all aircraft entering and departing a terminal area (where the greatest icing threat 
exists due to altitude and aircraft configuration constraints) and equipping some aircraft with an airborne 
capability to cover regions between terminal areas is seen as the best possible solution. Conceptual 
drawings of the ground-based and airborne capabilities are shown in Figures 1 and 2. 

The work described in this paper is being pursued under the NASA Aviation Safety Program’s 
Integrated Intelligent Flight Deck (IIFD) Project. Specifically, the research and development of remote 
detection and measurement of aircraft icing hazards is a part of the IIFD External Hazards Detection 
(EHD) element. While we are showing conceptual user interfaces here, the scope of the work being 
reported is focused on the sensing and characterizing of the external icing environment. EHD and IIFD 
have defined tasks to take the output of such sensing systems and determine expected flight path hazards, 
provide guidance to flight crews, and share the information to the airspace system. Since these sensor 
systems should provide extensive real-time environment assessment data, further linkage to FAA and 
NWS forecasting and nowcasting systems is also anticipated. 



Figure 1. — Ground-based remote sensing concept. 
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Figure 2. — Airborne remote sensing concept. 


Icing 
Nonicing cloud 


Airborne 


2.0 Ground-Based Remote Sensing Development 

2.1 Background 

The development of a ground-based hazard detection algorithm testbed, the NASA Icing Remote 
Sensing System (NIRSS), was begun in 1997 (Ref. 5). Ground-based detection was the first area of 
interest because it is not constrained by the space/power/weight limitations of flight systems. It was 
decided that initial efforts would be made with ground systems to prove basic icing remote sensing 
feasibility and then move to flight systems as technical maturity and programmatic constraints allowed. 
The objective of the ground-based system development is to provide accurate detection and warning of 
in-flight icing conditions in the near-airport environment with proven, cost effective, off-the-shelf 
meteorological sensing hardware. 


2.2 Instrumentation 

NIRSS utilizes a 3000 series multichannel radiometer (Ref. 6) built by Radiometrics Corporation. 
Ground-based passive microwave radiometers have the ability to profile temperature and water vapor in 
the atmosphere. They can also directly quantify total integrated liquid water between the antenna and the 
top of the atmosphere, and obtain low resolution profiles of liquid water; however this resolution is 
generally inadequate for our application. Weather radars have good range resolution of reflectivity to 
determine where clouds are and give some information on their composition, but do not measure water 
vapor, temperature or liquid water content. There is therefore synergism between these two sensing 
platforms, provided that they can simultaneously sample the same volume. Microwave radiometer 
beamwidths are typically several times wider than those of weather radars, but are sufficient for most 
meteorological applications, allowing the use of more compact and economical antennas. However, to 
match the sample volume of weather radars, much larger apertures and accurate pointing are required. 

Radiometric temperature profiling can be accomplished by profiling the downwelling spectrum on 
either side of the broad assemblage of oxygen transitions centered on 60 GHz and spanning from 
approximately 50 to 70 GHz. Oxygen emits proportional to its physical temperature, and the varying 
attenuation from the center of this feature out onto the wings yields range (altitude) information. Water 
vapor profiling can be accomplished by measuring the pressure broadening of the water vapor emission 
line at 22 or 183 GHz. Because the shape of these lines is dependent upon local pressure, information on 
the pressure altitude is contained in the emission profile, and density upon the intensity of the line. Cloud 
liquid water has no emission lines, but emits over a broad range of frequencies roughly as the square of 
the frequency. By observing widely disparate frequencies across the microwave spectrum, integrated 
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Figure 3. — NASA’s icing remote sensing site. Includes the 
larger trailer which houses the X-band radar (under the 
radome) and the instrument computers, the Radiometrics 
3000 series radiometer at the top left, the Vaisala ceilo- 
meter to the left, and the Metek Ka-band radar housed 
in the smaller trailer at the bottom of the photo. 

liquid water (ILW) measurements can be obtained. Including climatology as a priori information 
increases the skill of obtaining these measurements. The radiometer also has a vertically-pointing sensor 
which measures the infrared temperature (IRT) of the cloud base. 

A modified Honeywell WU-870 airborne X-band radar in a vertically-staring configuration is 
mounted atop the trailer which houses the system’s data fusion hardware (Fig. 3). X-band radar was 
chosen over longer wavelengths because the systems are generally cheaper and antennas are smaller. 
X-band also allows some sensitivity to cloud-sized water drops or ice crystals while not becoming 
significantly attenuated as are shorter wavelengths. A Vaisala laser ceilometer completes the current 
system configuration. 


2.3 Software Logic 

The instrument data are ingested into three separate computers within the NIRSS trailer for display 
and storage. The three raw data streams are pushed to a Linux machine where the fusion processing and 
hazard detection software reside. The processing begins by combining the radar and ceilometer data to 
determine where cloud layers exist. Ceilometer data are relied upon to determine cloud base heights in 
cases of light precipitation. If the radar reflectivity (REFL) is greater than 1 dBZ above the minimum 
detectable threshold for at least 200 m depth, a cloud layer is detected. Shallower clouds are ignored. The 
radar and ceilometer agree reasonably well during most of the examined case studies. Cloud boundaries 
are often quite noisy so a 5 min smoothing function is applied at no expense to accuracy or safety. 

Through years of in situ analysis, it has been found that the vertical distribution of liquid water in 
cloud layers cannot be determined accurately and reliably, even with the aid of external clues such as 
temperature range or cloud depth. Unless a full-blown cloud microphysics model is run with good 
initialization data including appropriate cloud active nuclei concentrations, one would still need to make 
assumptions about the liquid water profile given previous measurements of liquid water and ice versus 
water composition. NIRSS employs a simple fuzzy logic approach to estimate the liquid water profile, 
which can easily be modified in future upgrades as more case studies are collected and analyzed. For each 
cloud layer we have cloud top, base and depth. For the whole vertical column we have ILW, as well as 
temperature and REFL distributions with height. ILW for each cloud layer is assigned based on cloud 
depth and the coldest temperature in the layer. 
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The algorithm assumes four possible distributions of liquid within the cloud: uniform (same LWC for 
every level); wedge (from zero at the base to a maximum at the top); weighted by temperature (assume 
colder clouds with mixed-phase conditions have less water and more ice as the temperature is lower); and 
weighted by REFL 1/2 (an approximation (Ref. 7) for nonprecipitating clouds). The final calculated profile 
will be a combination of these four distributions scaled by weighting factors (aunifoim, a We d g e, a Te mp, and 
aREFL) for each, which are characteristics of the cloud that are derived from the data. These characteristics 
and weighting functions can differ for each cloud layer. Many combinations of weights are possible and 
they will depend on what we can learn about the cloud structure with the NIRSS measurements. Next, a 
determination about the cloud layer’s phase is made. If a cloud layer’s base temperature is greater than 
260 K and its REFL is less than 10 dBZ above the minimum detectable, the layer is deemed to be liquid 
phase. Otherwise, the layer is categorized as mixed phase. The layer’s LWC profile is then calculated 
with the weighting functions shown in Table 1, which have been derived from ongoing analysis of icing 
research aircraft flights. Having determined the heights of clouds and having the LWC and temperature 
profiles within them, a qualitative in-flight icing scheme (Ref. 8) is employed assuming a 15-pm drop 
size. The NIRSS display shows severity categories as designated in Table 2. 


TABLE 1.— WEIGHTING FUNCTIONS FOR FOUR METHODS OF 
DISTRIBUTING LIQUID WATER IN CLOUD LAYERS FOR 
LIQUID AND MIXED PHASE CLOUDS 


Weighting function 

Liquid cloud 

Mixed cloud 

^■Uniform 

0.1 

0.15 

awedse 

0.3 

0.25 

aTemp 

0.1 

0.5 

aREFL 

0.5 

0.1 


TABLE 2.— ICING SEVERITY CATEGORIES AND 
CORRESPONDING LWC VALUES SHOWN 
ON THE NIRSS ICING HAZARD DISPLAY 


LWC (gm“ 3 ) 

Icing severity 

<0.01 

No icing 

0.01 to 0.017 

Trace 

0.017 to 0.03 

Trace-light 

0.031 to 0.066 

Light 

0.067 to 0.012 

Light-moderate 

0.02 to 0.2 

Moderate 

0.20 to 0.37 

Moderate-heavy 

>0.37 

Heavy 


2.4 Recent Progress 

Many important improvements have been implemented within NIRSS in the last several years. The 
fusion software has been recoded from Lab VIEW to C++ language, and the display software has been 
recoded in the Java language. Realtime processing and analysis of data are now possible within the 
system trailer, and processed data and imagery are now automatically pushed to the web (http://icebox- 
esn.grc.nasa.gov/RSData/index.html). 

One recent enhancement to the NIRSS logic software is an algorithm which adjusts the radiometer’s 
temperature profile in cases of significant temperature inversions. Previous research has shown that 
radiometers often do poorly in detecting inversions (Ref. 9) due to their low vertical resolution, and also 
that inversions often are associated with the upper extent of significant super-cooled liquid water (SLW) 
(Ref. 10). If inversions are not accurately detected, ILW might be spread over a larger vertical range, thus 
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decreasing the detected severity. The new algorithm tests if the IRT is greater than 3° less than the 
radiometer derived temperature at cloud base. If so, the profile is rebuilt at the moist adiabatic lapse rate 
from the IRT through the cloud depth. An inversion with a fixed lapse rate is created starting at cloud top 
and then the rebuilt profile is gradually merged back to the radiometer’s profile by 2 km above the rebuilt 
inversion. Another enhancement is the quantitative use of the REFL field to warn on small-drop icing 
situations (Ref. 11). 

Work is ongoing to compare the icing severity output fields from NIRSS to NCAR’s Current Icing 
Product (CIP) and pilot reports (PIREPs). It is important to note that PIREPs are not considered highly 
accurate indicators of icing severity as they are prone to recording displacements in time and space and 
are also subjective measurements by an otherwise busy pilot. They are, however, the best measure of in 
situ icing hazard available short of research aircraft data. PIREP icing severities are scaled as 0 for null, 

1 to 8 for trace to heavy. 

An icing case where NIRSS was deployed at Cleveland’s Hopkins International Airport (KCLE) from 
0000 UTC on 26 February through 1600 UTC on 27 February, 2008 is shown in Figure 4. Starting just 
before 0000 UTC on the 26th, CIP detected light icing up to 18 Kft 1 which builds gradually down toward 
the surface by 0700 UTC on the 26th. A moderate PIREP, shown as a vertical stack of 5’s on the upper 
CIP severity plot, was reported around 0300 UTC from 10 to 14 Kft. NIRSS diagnosed trace to light icing 
in the same altitude ranges as CIP from before 0000 to 0500 UTC, with light icing shown at the time and 
altitude of the moderate icing PIREP. At 0500 UTC, it began snowing at KCLE and continued through 
1200 UTC on the 26th. The snow production and precipitation processes in clouds are well known to 
scorn most available liquid water through accretion of liquid to ice particles. NIRSS had an obvious gap 
in icing hazard profile from 0500 to 1300 UTC during the times of significant precipitation rate at KCLE, 
but CIP icing severity actually increased to moderate, with some altitudes showing heavy severities. This 
is due to a known problem with CIP at night (through 1300 UTC for KCLE) where the algorithm loses its 
visible GOES input and relies too heavily on GOES IR channels. From 1300 to 2400 UTC, precipitation 
rates decreased to intermittent trace amounts. CIP detects moderate severities generally from 4 to 18 Kft 
during the period. NIRSS showed periodic trace severities from 1 to 14 Kft during the period. 



Figure 4. — CIP (A) and NIRSS (B) derived in-flight icing severity over KCLE from 0000 UTC on February 26, 2008 to 1600 
UTC on February 27, 2008. Precipitation rates recorded at KCLE are shown in part C. 


'Note that feet, rather than meters, are used in NIRSS since this is the unit convention used in aviation. 
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Two null PIREPs occurred just after 1400 and 1500 UTC at 7 to 8 Kft, respectively. A trace PIREP 
occurred over KCLE between 8 and 9 Kft just after 1800 UTC, which was coincident with light severities 
by both CIP and NIRSS. NIRSS later missed a light to moderate PIREP recorded somewhere between 1 
and 20 Kft just after 23 UTC on the 26th. From 0000 to 1600 UTC on the 27th, both CIP and NIRSS 
indicated a gradual lowering in heights of the icing layer and diminishing severity through the profile. 
Several null PIREPs were recorded between 14 and 16 UTC between 4 and 8 Kft that corresponded to no 
icing in NIRSS and trace in CIP. 

Overall, NIRSS and CIP did well in diagnosing similar severities to PIREPs for the case. Both 
products saw the case evolve similarly in time and altitude ranges. NIRSS seemed to under predict the 
real extent of hazard in time, but did well with the hazard altitude ranges. NIRSS also seemed to 
somewhat under predict the severity value when ILW was detected. Both of these deficiencies could be 
easily scaled if further case studies confirm the biases. CIP seemed to exhibit an over prediction in 
severity value at any given time, especially at night when the most significant snow scavenging was 
known to be occurring. Both the NIRSS and CIP products should be able to be improved upon by case 
study inter-comparisons such as the one in this work. 

2.5 Future Work 


2.5.1 Software Upgrades 

Several possible upgrades to the existing NIRSS logic will be explored in the future. Examining the 
three (radar, ceilometer, radiometer) cloud base estimates under different weather conditions may yield 
additional information on precipitation intensity or type that we can use further within the system. 
Additional analyses of aircraft research flight data along with NIRSS input data should lead to 
improvements on the existing simple algorithm for ILW distribution. Another improvement that is needed 
is refactoring of the C++ codebase to allow for the parameterization of keyword inputs and allow future 
system hardware generalization. The NIRSS logic code is a body of work in progress, and many 
improvements have been implemented in the recent past. A software redesign would make the package 
more flexible for recent and future enhancements. 

NIRSS ’s severity definition is based on the amount of time it takes an aircraft to accrete 1/4 in. of ice 
on the leading edge from less than 5 min which is classified as “severe”, to greater than 15 min, which is 
classified as “light”. Work needs to be done to see if adjustments are needed to bring aircraft icing time 
scales into accord with the severity scales that are being measured with NIRSS. More information could 
be extracted from the system to make a better estimate of drop size, or at least flag when SLD conditions 
are expected. This could be accomplished by using vertical profiles of radar REFL along with the rain 
flag from the radiometer to detect bright bands. 

Another important software upgrade would be the inclusion of a module that ingests the closest CIP 
gridpoint data and NWS balloon sounding data. The possible synergies between the ground-based NIRSS 
and the model/satellite based CIP icing hazard products should be fully utilized by NIRSS. 

2.5.2 Instrumentation Upgrades 

One area where significant improvement is possible is in determining cloud phase. The Radiometrics 
89/150 GHz radiometer could help determine when a cloud is mixed phase by examining the ratio of 
reflected brightness temperatures at the two frequencies. Previous research (Ref. 12) with this channel 
combination has shown that ice water path and effective particle diameter can be calculated at any 
viewing angle. Radiometrics Corporation is currently building a scanning “pencil-beam” radiometer for 
NASA Glenn with channels from 22 to 30 GHz and at 89 GHz. The inclusion of vertical and horizontal 
polarizations at 89 GHz allows for particle phase discrimination as SLW droplets are axial-symmetric and 
therefore have no preferential scattering polarization. Glaciated particles tend to grow and fall 
preferentially in one axial dimension, thus creating a larger difference in the polarization channels. 
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Upgrading the radar to a Doppler X-band would provide information on mixed phase conditions by 
being able to classify fall velocities spectra as rimed ice particles, unrimed ice particles or small 
suspended particle populations (Ref. 13). Rimed particles unequivocally indicate the presence of 
supercooled liquid. 

2.5.3 Research Flights 

There are good data sets from the NASA Twin Otter, the MSC Convair and the UND Citation from 
the Second Alliance Icing Research Study (AIRS-II) field project in 2003 and Cleveland-area flights 
during 2005. Additional research flights with the new NASA S-3 would be extremely valuable to 
facilitate future NIRSS improvements. Now that the system has been proven to accurately detect icing 
hazard in its vertically staring configuration, trials should begin to compare the scanning pencil-beam 
radiometer and the 3000 series multichannel radiometer to the nearby KCLE NEXRAD output and any 
future icing research aircraft flights. The power of NIRSS to warn of icing hazard in the near airport 
environment will be greatly enhanced when elevation and azimuth scanning capabilities are implemented 
and tested. 


3.0 Airborne Remote Sensing Development 

3.1 Airborne Radar Development 


3.1.1 Background 

Recent studies (Refs. 14 to 16) concluded that colocated microwave and millimeter wavelength radars 
can estimate liquid water content and a characteristic drop size of liquid precipitation, suggesting the 
application of the technology for in-flight sensing of icing potential ahead of aircraft. Although there are 
some differences between the various retrieval techniques, all the algorithms estimate moments of the 
drop size distribution from profiles of reflectivity at a weakly attenuating microwave frequency, typically 
X-band, and a more attenuating microwave frequency, 35 and/or 95 GHz. The microwave channel is 
highly sensitive to drop size, as the reflected power is proportional to the sixth moment of the drop size 
distribution (p r = f(IDj 6 )) (for Rayleigh scatterers), while the differential attenuation, between the 
microwave and millimeter wave channels, can be related to liquid water content. After successful ground- 
based tests, the Icing Branch of the NASA Glenn Research Center began collaborating with ProSensing 
Inc. to develop an airborne, multifrequency radar to test the technique in realistic in-flight icing 
conditions. 

3.1.2 Instrumentation 

The NASA Multifrequency radar (MFR) system was designed with 2 mm wavelength channels 
(35 and 95 GHz) to extend the liquid water range of the sensor and to maintain accuracy in nonRayleigh 
precipitation. The radar consists of pulsed X (3 cm) and Ka-band (1 cm) radars, both employing 25 kW 
magnetron transmitters, and a two-antenna Frequency Modulated Continuous-wave (FMCW) W-band 
(3 mm) radar. With the RF sections housed in a single forward looking pod, the three radars will sample 
similar sample volumes with 12, 5.5, and 3 in. diameter antennas at 10, 35, and 95 GHz, respectively. The 
95 GHz antennas are Gaussian optics lenses, while the Ka and X-band antennas are flat-plate slotted 
waveguide arrays (left side of Figure 5). When installed on the NASA S-3 aircraft, all three forward- 
looking radars will transmit and receive horizontally polarized signals. The front section of the pod, 
shown on the right of Figure 5, contains the antennas and a hemispherical X-band radome and the Ka and 
W-band Rexolite radome windows tilted at the Brewster angle for maximum signal transmission. 
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W-band 

3 in. lens antenna 


| 


W and Ka-band Rexolite windows, 
tilted at the Brewster angle 
(58° at H-pol) 



Ka-band 5.5 in. 
slotted waveguide 


X-band 12 in. slotted waveguide antenna \ 

Figure 5. — The NASA Multifrequency Radar (MFR) system pod section, containing an X (9.4 GHz) and 
Ka-band (35 GHz) pulsed radars and a W-band (95 GHz) FMCW radar system. 


3.1.3 Data Analysis 

In a 1998 study (Ref. 17), a working hypothesis that cloud parameters could be extracted by 
measuring backscatter at a combination of attenuating and nonattenuating frequencies was developed. 
Since scattering is a complex nonlinear function of particle size and frequency, it is impractical to 
consider an analytical solution to the inverse problem of computing particle size and liquid water content 
based on measured backscattered power at multiple frequencies. The analysis effort for the MFR is 
therefore focused on an approximate numerical solution to the inversion, specifically, a neural network. 
Neural nets are ideal for solving problems where the forward problem is well characterized but the 
inverse is nonlinear and complicated. The neural net programming process starts with the generation of a 
training set of statistically significant atmospheric conditions. This took the form of profiles of drop size 
distributions. From these drop distributions, a computer algorithm calculated the corresponding radar 
observed reflectivity profiles at each operating frequency as well as resultant LWC and drop size in each 
volume cell. Gaussian distributed random variables were added to the reflectivity profiles to simulate the 
effect of noise in the data. This set of multifrequency reflectivity profiles (inputs) and cloud and 
precipitation parameter profiles (outputs) were then used to train the neural network (Ref. 15). The 
network was trained by simulating thousands of test cases of radar scattering from assumed particle size 
distributions. 

3.1.4 Progress 

X, Ka, and W-band data from the Mt. Washington Icing Sensors Project (MWISP) (Ref. 18) were 
analyzed with the neural net. Figure 6 shows the comparison of the neural net retrieval’s LWC profile to that 
measured with the ATEK LWC profiling radiosonde. The astericks show the average LWC as derived by 
the neural network. These points were computed by taking 12-min averages of the neural network derived 
LWC at each altitude. Figure 6 shows good agreement between the in-situ radiosonde derived LWC and the 
neural networks output. The measurements show agreement both in the altitude distribution (within the 
resolution of the measurement) and peak value (within 20 percent). Further time history plots and the 
corresponding weather observations are included in Pazmany’s 2001 report (Ref. 15). 


NASA/TM— 2009-21 5828 


9 



Liquid water, g/m 3 

Figure 6. — Coincident 10, 35, and 95 GHz Radar derived liquid 
water profile from MWISP field project marked with *. Comparison 
to ATEK radiosonde in-situ LWC measurements marked with line. 


3.1.5 Future Work 

The MFR system is expected to be delivered to NASA following modification in the latter part of 
2008. Following delivery to NASA, testing is planned in the NASA Icing Research Tunnel (IRT) to 
assess the ice accretion patterns on the X-band radome and the Ka and W-band windows and to assess its 
general airworthiness. Following this IRT testing and vertical staring ground-based evaluation of the radar 
system, the MFR will be installed and tested on the NASA S-3 aircraft. 

3.2 Microwave Radiometry 


3.2.1 Concepts 

A passive airborne approach has been developed by the Army Corps of Engineers Cold Regions 
Research and Engineering Laboratory (CRREL) using microwave radiometers to detect icing conditions 
ahead of aircraft. Millimeter wave (MMW) radiometry offers the potential of smaller size, less weight, 
lower power consumption and lower cost than active systems. CRREL is working with systems operating 
near 35 and 94 GHz. 

Microwave radiometers operate by receiving thermal energy from atmospheric constituents at specific 
wavelengths in the electromagnetic spectrum. Passive microwave radiometers traditionally measure 
atmospheric and terrestrial characteristics from ground-based, such as NIRSS, and satellite-borne 
systems. However, their ability to measure atmospheric temperature, cloud liquid water content, and 
attributes of cloud and precipitation constituents, especially viewing along a near-horizontal flight path, is 
less well developed. The radiation intensity observed by a radiometer is a function of the temperature, 
reflectivity, transmissivity, and emissivity of the emitter and attenuation by constituents between the 
emitter and the radiometer at the wavelength of interest. 

Gas, liquid, and solid atmospheric constituents each have unique absorption and emission spectra due 
to their dielectric properties. As a result, the atmosphere absorbs in several narrow wavelength bands and 
allows radiation to be transmitted through several broad windows. The primary absorbers are oxygen and 
water vapor. Oxygen absorbs and reemits in the 50 to 60-GHz region and at 1 18 GHz and is used for 
temperature profiling. Water vapor has peak absorption and reemission at 22, 37, and 183 GHz (Ref. 19). 
Peak absorption and emission of liquid water occurs near 37 and 89 GHz. Though the dielectric strength 
of water is sufficiently high that is absorbs and emits strongly, the frequency-dependent dielectric strength 
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of ice is approximately 10 percent of that of water in the millimeter wave spectrum. Therefore, to a 
MMW radiometer the brightness temperature of water clouds of a given water content tend to be greater 
than the brightness temperature of ice clouds of a similar amount of frozen water. Therefore, supercooled 
water clouds are generally more easily detected by microwave radiometers, an advantage because they 
present the greatest hazard. Ice crystals typically do not adhere to aircraft wings and therefore do not 
contribute to the icing hazard. However, because water vapor absorbs and emits near 37 GHz and is 
potentially distributed nonuniformly in the clear atmosphere outside the cloud, the behavior of sensors 
will also vary within and between flight levels (Ref. 20). 

In addition to the absorption, emission, and scattering of atmospheric constituents, the Earth’s surface 
and space both have different signatures in the MMW spectrum. Depending upon materials at the Earth’s 
surface, such as land or water, the Earth’s surface will typically provide a large brightness temperature at 
35 and 94 GHz because the surface is warm. Conversely, the cold of space provides a small brightness 
temperature at 35 and 94 GHz. These characteristics may be used to advantage to detect conditions ahead 
of an aircraft through contrast that they provide with the brightness temperature of clouds. 

Savage (Ref. 20) and CRREL developed a concept that will allow cloud liquid water content to be 
derived for clouds, and may provide an estimate of cloud temperature from brightness temperatures. 
Radiometers placed on the nose of an aircraft could scan ahead of the aircraft horizontally, and 2° above 
and below horizontal. In a clear-sky condition, the 2° beam detects colder temperatures by viewing cold 
space than does the -2° beam looking toward the warmer surface of the Earth. Since the atmosphere is 
more transmissive at 35 GHz, the brightness temperature of the 2° beam will be less than (colder) than the 
2° 94 GHz beam. For the downward looking geometry, the brightness temperature difference between 
35 and 94 GHz depends on the surface and the atmospheric conditions and is less easily interpreted. The 
horizontal beam may provide an indication of cloud sensible temperature. 

In operation, as the aircraft approaches a cloud, the 2° beam and the -2° beam converge in brightness 
temperature toward that of the horizontal beam. The rate of convergence to the temperature of the 
horizontal viewing beam depends on the distance from the cloud, the cloud height above and below the 
flight altitude, the decrease of air temperature with altitude, and the cloud liquid water content. An 
estimate of liquid-water content magnitude is obtained by comparing the brightness temperatures of the 
35- and 94-GHz beams in the 2° orientation. Savage et al. (1999) also speculated that the presence of 
drizzle-size drops can be detected by sensing polarized radiation resulting from drizzle drop scattering of 
polarized radiation from Earth surfaces, and shape distortion of the largest drops. The 2° upward-looking 
and horizontal-looking radiometer concepts may provide the most useful information with regard to 
locating supercooled water because it does not involve the complexity of changing brightness 
temperatures emitted or scattered from the Earth’s surface. 

3.2.2 Simulation 

We modified the RADTRAN model, developed by the former Air Force Geophysics Laboratory, to 
model attenuation, scattering and brightness temperatures for 1 to 1000 GHz (Ref. 21). Originally 
developed to predict brightness temperature and attenuation for ground based and space-based microwave 
radiometers looking zenith, nadir, or at an angle, RADTRAN calculates the attenuation, transmittance, 
and brightness temperatures for thermal emission at microwave frequencies using six predefined 
atmospheric profiles. A capability to model partial or fully glaciated clouds is not available in 
RADTRAN, though precipitation can be either in a liquid or solid state. 

A virtual fly-through capability was created in RADTRAN for simulating in-flight icing. The user 
specifies the cloud and atmospheric conditions, the frequencies and view directions, the start and end 
point of the flight path, the number of points along the flight, and the altitude of the flight path. At 
designated points along the flight path the cloud atmospheric conditions along the view direction are 
specified. RADTRAN then computes the brightness temperatures, attenuation, and transmission for each 
frequency and view direction along the flight path location. Figure 7 presents an example of a 37 GHz 
radiometer with 90° and 92° viewing angles used to compute both the distance to a cloud ahead of the 
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Figure 7. — Brightness temperatures associated with 90° and 92° view 
elevations for radiometers on an aircraft approaching a single cloud. 


aircraft and the cloud top altitude for a single cloud with a uniform cloud liquid water content. From the 
changes in brightness temperature for the two viewing angles, aircraft altitude and the arc distance S 
between when the 90° viewing sensor first detects the cloud and the 92° viewing sensor detects the cloud, 
it is possible to compute the cloud top altitude and the distance to the cloud from the position when the 
92° viewing sensor first detects the cloud. The calculated cloud top is 0.96 km and the distance to the 
cloud is 9.2 km, while the actual values are 1 .0 and 1 0 km, respectively. 

Massie et al. (Ref. 22) have developed a neural network model to predict liquid water content for the 
CRREL system. The net was trained using approximately 10,000 cases and evaluated using independent 
cloud scenarios. The model output consists of liquid water content and brightness temperature as a 
function of frequency and view direction. Predicted liquid water contents agreed very well with the 
observed liquid water contents. 

3.2.3 Hardware Development 

Two flight-ready polarimeters developed by WaveBand and Sierra Nevada, Inc. are now available for 
flight testing of the icing condition detection concepts modeled in RADTRAN and with the neural 
network. The polarimeters provide the four Stokes parameters enabling the potential detection of drizzle 
size and larger drops, and ice crystals. Initial flight tests have been conducted but explicit corroborative 
cloud microphysical information was not available for those first flights. These or similar frequency 
sensors are also currently planned for flight on the fully instrumented NASA-GRC icing research aircraft 
within the next several years. 


4.0 Remote Sensing Validation Method Development 

4.1 Background 

Verification and validation of the accuracy of the remote sensing instrumentation requires concurrent 
in-situ measurement of the atmospheric conditions of interest. A large suite of suitable instruments have 
been developed for use on aircraft and provide valuable V&V data. However, access to instrumented 
research aircraft is limited and another method for in-situ characterization of the atmosphere has been 
identified. Standard, commercially available balloon-borne radiosondes are used to provide a limited 
number of atmospheric parameters, specifically temperature and relative humidity as a function of altitude 
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and position. These sondes lack the ability to measure liquid water content, which is vital to determining 
the level of severity of the icing hazard. NASA has undertaken an effort to develop a LWC sensor which 
is suitable for use on disposable radiosonde. 

4.2 Instrumentation 

The radiosonde system used by NASA is a Sippican Model W-9000 Meteorological System, an 
integrated hardware and software package that provides meteorological profiles through the use of 
telemetered radiosonde data and GPS location data. Temperature and humidity are directly measured, 
while pressure is calculated from GPS data and the Standard Atmosphere. Temperature, with an accuracy 
of 0.2 °C is obtained using a standard rod thermistor, positioned at the end of an arm to remove it from 
the close proximity of the sonde case. A deposited carbon film on an acrylic substrate is used for humidity 
measurement with a response time of typically less than one second. Comparison of GPS positional data 
from the sonde and the ground station is used to calculate sonde altitude, latitude and longitude, its 
velocity with respect to the launch site and rate of ascent. We have recently upgraded this system to the 
Lockheed Martin Sippican Model LMG6 which uses the same sensor technology but with improved data 
acquisition signal processing. Data from radiosondes were used during the MWISP, AIRS I, and AIRS II 
field studies to develop the remote sensing instruments. Figure 8 shows a comparison of data between the 
Radiometrics 3000 radiometer (blue trace) and the radiosonde (red trace) taken during AIRS II (Ref. 23). 
Temperature agreement is very good; that of relative humidity is good, with the exception of those 
regions where the RH is rapidly changing. However, the radiosonde cannot measure nor provide an 
estimate of LWC, the right-most plot in Figure 8. 


4.3 Future Work 

Several approaches have been used to develop a liquid water content probe of sufficiently low mass, 
power consumption and cost for use on a disposable radiosonde. The ATEK device, which uses a 
vibrating wire to measure LWC and requires temperatures below freezing, has been flown with some 
success. However it is no longer available as it never reached commercial status. Innovative Dynamics, 
Inc., through a SBIR Phase II contract, is currently developing a probe based on optical sensing. Two 
near-IR laser sources operating at wavelengths of 1310 and 1550 nm are used to illuminate the cloud. 
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MP: 18:14:34 11/11/03 TIR: 271.8 K 


Rain: Y 



10 

9 

8 

7 

6 

5 

4 

3 

2 

1 




Int Vap: 1.91 cm 
Int Liq: 0.32 mm 


10 

9 

8 

7 

6 

5 

4 

3 

2 

1 


Relative humidity Liquid, g/m 3 


Figure 8. — Comparison of remotely sensed and in-situ data from AIRS II field study. 
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The intensity and polarity of the backscattered light is used to determine LWC, droplet size and density 
and to discriminate between liquid and ice. The concept has been demonstrated in breadboard form and 
has been extensively tested in a cloud chamber. The prototype flight article is nearing completion, initial 
testing will be in a tethered balloon configuration. The use of a small, low speed, vertical icing tunnel 
under development at Glenn is planned to calibrate the probe. Deliverables include three flight-ready 
probes that will be flown using the Glenn radiosonde system. 
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